Slow Logarithmic Decay of Magnetization in the Zero Temperature Dynamics of an 
Ising Spin Chain: Analogy to Granular Compaction 
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We study the zero temperature coarsening dynamics in an Ising chain in presence of a dynamically 
induced field that favors locally the ' — ' phase compared to the '+' phase. At late times, while the 
'+' domains still coarsen as t^^^, the ' — ' domains coarsen slightly faster as t^^^ \og{t). As a result, 
at late times, the magnetization decays slowly as, m{t) = —1 + const. /log(t). We establish this 
behavior both analytically within an independent interval approximation (IIA) and numerically. 
In the zero volume fraction limit of the '+' phase, we argue that the IIA becomes asymptotically 
exact. Our model can be alternately viewed as a simple Ising model for granular compaction. At 
late times in our model, the system decays into a fully compact state (where all spins are ' — ') in a 
slow logarithmic manner ~ l/log(f), a fact that has been observed in recent experiments on granular 
systems. 
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The effect of quenched disorder on the relaxation dy- 
namics of many body systems has been studied quite ex- 
tensively In systems such as structural glasses, where 
quenched disorder is absent, an alternative approach has 
been put forward that considers the slow relaxation due 
to kinetic disorder, induced by the dynamics itself 
Another important system where kinetic disorders give 
rise to slow relaxation is granular material. The density 
relaxation of loosely packed glass beads has been studied 
in recent experiments and it was found that the density 
p{t) compactified slowly as, p(oo) — p{t) ^ l/log(f), un- 
der mechanical tapping It is natural to expect that 
such kinetic disorders may play important roles in the dy- 
namics of other systems as well. In this Letter we study, 
for the first time, the effect of a dynamically self-induced 
field in an important class of out of equilibrium problems, 
namely the domain growth problems, and show that such 
systems also exhibit logarithmic relaxation, suggesting 
that inverse logarithmic relaxation is quite robust. 

Domain growth following a rapid quench in tempera- 
ture in ferromagnetic spin systems is one of the better 
understood out of equilibrium phenomena . For exam- 
ple, if an Ising system is quenched rapidly from infinite 
temperature to zero temperature without breaking the 
symmetry between the two ground states, domains of up 
and down spins form and grow with time. The aver- 
age linear size of a domain grows with time as l(t) ^ 
for zero temperature nonconserved dissipative dynamics. 
However if one puts on a small uniform external field (say 
in the down direction) , then even at zero temperature the 
system quickly reaches the pure state of magnetization 
— 1 in a finite time proportional to the initial size of the 
up domains. A natural question is: what happens when, 
instead of applying a global external bias, the symmetry 
between the pure states is broken locally by the dynamics 



itself ? 

In this Letter we address this question in the con- 
text of a simple Ising spin chain with spins Si =^ ±1. 
Starting from a given initial configuration, the system 
evolves by single spin fiip continuous time dynamics. 
Let W{Si; Si-i, Si+i) denote the rate at which the flip 
Si —> —Si occurs when the two neighboring spins are 
Si-i and S'i+i. In our model the rates are specified as 
follows: 



-) = W{+;- 

-) = 1 
-) ~ a 



-) = o 



VF(- 



_ 1 

" 2 
(1) 



Here we restrict ourselves to the case a = 0. We note 
that the case a — 1 corresponds to the usual zero tem- 
perature Glauber dynamics The only difference is 
that for a = 0, the move (-1-,—,+) {+,+,+) is not 
allowed and thereby the symmetry between and '— ' 
spins is locally dynamically broken. Thus isolated '— ' 
spins (surrounded on both sides by a '-I-') block the co- 
alescence of domains and locally favor the '— ' spins. 
As a result, we show below, the system eventually decays 
into the state where all spins are '— ' but does so in a very 
slow manner, m{t) + 1 ^ l/log(i). 

Our main results can be summarized as follows. In con- 
trast to the case a — 1 (where the average size of both 
and '— ' domains grow as l±{t) ^ t^^^ at late times [Q and 
the average magnetization m(t) = {1+ — l_)/{l^ + 1_) is 
a constant of motion 0), for a = we show that at late 
times, whfle l+{t) ^ V^'^ , ^ t^^'^ log(6t) where 6 is a 
constant that depends on the initial condition. Thus due 
to the dynamically generated local bias, the '— ' domains 
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grow slightly faster than the '+' domains and as a result 
the magnetization decays as, m(t) = — 1 + const/log(6i) 
for large t. Notice that the average domain size grows 
faster for a = than for a = 1, i.e. paradoxically coars- 
ening is enhanced by putting one of the rates to zero. 

Our model can alternately be viewed as a toy model 
of granular compaction if one identifies the '— ' spins 
as particles,'-!-' spins as holes and the 1-d lattice as a 
section of the bottom layer of a granular system. The 
final state where all spins are negative (magnetization, 
m = — 1) then corresponds to the fully compact state 
with particle density 1. The basic mechanism for gran- 
ular compaction can be summarized as follows. There 
exist local kinetic 'defects' and the system can gain in 
compaction only by relaxing such local defects. Such 
relaxation happens via the tapping process. However, 
these defects become rarer with time and it becomes 
harder and harder for the system to find such a local 
defect, relax it and thereby gain in compaction. This 
is the origin of the slow logarithmic relaxation. In our 

model, the triplets 'H h' play the role of such local 

defects which decay only via the diffusion of kinks. Thus 
the diffusion effectively plays the role of tapping. The 
density of these triplets decays with time and the system 
finds it progressively harder to relax. Such logarithmic 
relaxation has been found previously in various models 
of granular systems Q|. However our model differs from 
these models in an important way. In previous models, 
the density, while increasing as p{oo) — p{t) ~ l/log(i) at 
intermediate times, finally saturates to its steady state 
value exponentially fast Q for finite values of the rates 
of microscopic processes. In contrast, in our model, the 
approach to the final state is logarithmic asymptotically 
even at very late times. 

In terms of the motion of the domain walls between 
'-I-' and ' — ' phases, our model can also be viewed as a re- 
action diffusion process. We need to distinguish between 
the two types of domain walls —+ = A and H — = B. 
Note that by definition (originating from a spin config- 
uration) the A's and i3's always occur alternately. The 
A's and B's diffuse and when an A and a B meet, they 
annihilate only if A is to the left of i?, otherwise there is 
hard core repulsion between them. 

To start with, we set up our notations. We define 
Pn{t) and Rn{t) to be the number of domains of size n 
per unit length of '-f ' and '— ' types respectively. Then, 
N{t) — J2„ Pn — X^n number of domains of 

either '-(-' or '— ' spins per unit length. The density of 
kinks is therefore 2N{t). We also define the normal- 
ized variables, Pn = Pn/N and r„ ~ Rn/N. pn (or 
r„) denotes the conditional probability that given a do- 
main of '-I-' (or '— ') has occurred, it is of length n. Let 
L+{t) — ^ nPn and L^{t) = ^ ni?„ denote the densities 
of '-I-' and '— ' spins. Clearly L^{t) + L^{t) ~ 1 and the 
magnetization per unit length is m{t) = L^{t) — L^{t). 



The average size of a and a '— ' domain are de- 
noted respectively by l+{t) = 'Y^npn = L+(t)/N and 

Following Glauber's calculation for the a = 1 case, it is 
easy to show that for a = case, the domain density 
N{t) = (1 — (S'iS'i+i))/4 of either phase, and the fraction 
of '-I-' spins, i+(i) = (l-f {Si))/2 evolve according to the 
exact equations. 
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where Pi{t) = ((1 - 5,_i)(l + S,){1 - S,+i))/8 and 
i?i(i) = ((l + 5,_i)(l-S',)(l-h5'j+i))/8 are respectively 
the density of '-I-' and '— ' domains of unit length, i.e., 

the density of triplets ' \ — ' and 'H h. It is easy 

to see physically the origin of these two exact equations. 
Eq.(||) arises from the fact that the domain density can 
decrease only via the annihilation of the triplets ' — I — '. 
Also, on average, the fraction of '+' spins can decrease 

only due to the blockage by '-I h' triplets giving rise 

to Eq.(|). 

Using m{t) = {S^) = 2L+{t) - 1, we find from Eq.(|) 
that dm/dt = — 2i?i. We note that for the case a — 1, 
dm/dt = indicating that the magnetization does 
not evolve with time. In our case, due to the triplet 

defects 'H h', the average magnetization decays with 

time. We also note that, unlike the a — 1 case, the evo- 
lution equation (^ for the single point correlation func- 
tion involves two and three point correlations (via Ri{t)). 
Writing down the analogous equations for i?„ (t) gives an 
infinite hierarchy which makes an exact solution difficult 
for a = 0. 

Using Ri — riN and = l^N in Eq.(||), one can 
formally solve for N{t) in terms of ri and l-^ as. 
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Furthermore if the density of the '-I-' phase is L^{tf)) ~ e, 
then, using the relation TV (t) = l/[l_{t)+l+{t)] inEq.(|), 
we find 



l-(t) 



■ exp 



' nit') 
. hit') 



dt' 



(5) 



clearly showing that the ratio l^{t) /l^{t) is growing due 

to the presence of the triplets 'H h'. Note that the 

asymmetry between the growth of ' — ' and '-I-' domains 

is evident due to the triplet defects '-] (-', present with 

density Ri — riN. 

In order to make further analytic progress, we first use 
the IIA where correlations between neighboring domains 
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are neglected. The IIA was used previously for the pure 
Glauber-Ising model, i.e., the a = 1 case It yielded 
results in agreement, qualitatively as well as quantita- 
tively to a fair degree of accuracy, with the exact results 
available [^j9| . Following the derivation of the IIA equa- 
tions in the a = I case, it is straightforward to derive 
the corresponding equations for the a — case . Un- 
der this approximation, the domain densities Pn{t) and 



Rn{t) evolve as |10 
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dt 
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for all n > 1 with Pq = (absorbing boundary condition) 
and 

dRn Pi 

~ Rn+l + Rn-1 — '2Rn — ^(^" + ^n-l) 



Pl_ 



n-2 



^R,Rn-^-l\ n>2 
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where N{t) = 'Y^Pn = X^^n- '^^'^ tie easily checked 
that these two IIA equations satisfy Eqs.(p|) and (||) ex- 
actly, and consequently also Eqs.(^ and (g). 

To calculate N{t) using Eq.(||), we need to evaluate 
two quantities from the IIA equations: (i) ri{t) = Ri/N 
and (ii) l+(t) — ^np„ where p„ = Pn/N. In order 
to calculate these two quantities, it is useful to write 
the IIA equations in terms of the normalized variables, 
Pn = Pn/N and r„ = i?„/iV. From Eqs.(|) and (0), we 
then get, 



dpn 
dt 



^ Pn+l + Pn-l - '^Pn + f 1 (Pn " Pn-l) + PlPn (8) 



for all n > 1 with po ~ (absorbing boundary condition) 
and 

—7- = rn+i + r„_i - 2r„ - _pir„_i 
dt 

n-2 



+Pi'^rirn-i-i; n>2 
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= r2- n. 



(9) 



It is easy to check that the normalization condition 
^Pn = X] ~ 1 satisfied by these two equations. 

The two IIA equations above are coupled nonlinear 
equations with infinite number of variables and hence 
are difficult to solve exactly. Our approach is a combina- 
tion of a scaling assumption and then rechecking this as- 
sumption for self-consistency. Consider first the p„ equa- 
tion, i.e. Eq.(|). We substitute p„(t) = t-^^^ f{nt^^/^ ,t) 
in Eq.(^) and ask if the resulting equation allows for a 
steady state scaling solution as t — > oo, i.e., if the scaling 



function becomes explicitly independent of t as < — > cx). 
It is easy to verify that if ri{t) decays faster than t~^^^, 
such a time-independent scaling solution is possible with 
f{x) — I exp(— a;^/4). In this case, l+{t) = "Pn ~ 



t^^^ /o°° xf{x)dx = \fwi at late times. 

Next we consider the r„ equation, i.e., Eq.(H). Since 
P\ = —dlogN/dt, a natural choice would be to write 
r„(t) = N{t)g{nN(t),t). Substituting this in Eq.(|), we 
find that in the t —* oo limit, the equation allows for a 
time independent scaling function, g{x) = ccxp(— ca;) (c 
is a constant), provided N{t) decays faster than In 
this case, ri = N{t)g{0) = cN{t). Thus if scahng starts 
holding beyond some time to, then c = ri{tQ)/N{to). 

Using the results (i) ri{t) = ri{to)NU)/N{to) and (ii) 
1+ {t) = spivt in the exact equation Eq. (Q) , we find 



where log(5) 



N{to) 



log(&) 
t \ogibt/to) 



(10) 



ri(to)\/to 

expression for ri (t) , we find 



riit) = 



Substituting this result in the 



Vtlog{bt/to)' 



(11) 



We now use the late time result Eq.(|T^) in the exact re- 
lation Eq.(0) and find, 



Pi 



2t tlog{bt/ta)' 



(12) 



From the above expressions, it is clear that both ri(t) 
and N{t) decay faster than t~^/^ and hence our scaling 
solutions are completely self-consistent. 

It is easy to see that these IIA results become exact in 
the zero density limit of the phase (e ^ 0). In this 
limit, the average size of a ' domain is 1/e times larger 
than the average size of a domain. As time increases, 
the domains will certainly grow in size. But a typical 
'-f ' domain will disappear (via the absorbing boundary 
condition) much before encountering other '-f ' domains, 
i.e., before feeling the presence of the constraint due to 

triplets '-I h'. The probability of such an event is of 

order 0(e). Thus, effectively, the dynamics of the system 
will proceed via eating up of the domains. Hence, 
if there is no correlation between domains in the initial 
condition, the dynamics is not going to generate correla- 
tions between them and hence IIA becomes exact. The 
picture in this limit is similar to the zero temperature 
dynamics of the q state Potts model in the limit g ^ 1+ 
Q. For other volume fractions, it is likely that IIA pre- 
dicts the correct fixed point picture at late times. This 
is confirmed by Monte Carlo simulations of the model. 

To improve efficiency, these simulations were made for 
a version of the model with simultaneous updating. This 
should not change any of the above conclusions. For 



3 



convenience, we chose initial conditions such that all do- 
mains of minority spins had length 1, while the lengths 
of majority spin domains were distributed exponentially. 
At each time t, all kink positions were written into an 
array K, and only this array is used to generate the array 
K ' for the next time step. For each value of m(0) we sim- 
ulated between 20 and 200 lattices of 2^^ sites for 3 x lO'' 
time steps. 



ditions is the current value of the magnetization). But 
the detailed predictions for the scaling function show sub- 
stantial corrections which seem however to disappear for 
t oo. 
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I (0) / 1.(0) = 20 
1(0)/ 1,(0) = 14 
1(0)/ 1.(0) = 10 
I (0) / 1, (0) = 5 
I (0) / 1, (0) = 2 
L(0) / 1 (0) = 1 
I (0) / 1 ,{0) = 0.3 
1(0)/ 1,(0) = 0.1 
l{0)/U{0) = 0.03 
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FIG. 1. l-{t) /l+{t) versus t for runs with different initial 
magnetizations. 



FIG. 2. y/tN{t) versus t/T[m(0)] for the same runs as in 
Fig.l. The solid line is the predition const/ log(t/r). 

Data for l-{t)/l+{t), plotted in Fig.l, show the pre- 
dicted monotonic increase with t. This increase is log- 
arithmic for t > to, while it is faster for t < to. For a 
detailed comparison with the above theory we need to 
know how ta (and thus also b) depends on to(0). We 
expect it to be exponential for m(0) << 1, but this is 
not sufficient for a detailed analysis. Thus we determine 
for each m(0) a t such that the data for for 
^/tN{t), for tpi{t) and for ^/tri{t) collapse when plotted 
against t/r. The fact that a single T[m(0)] exists which 
gives a good data collapse in all four plots is highly non- 
trivial. We just show such plots for N{t),pi{t) and ri(t) 
in Figs. 2 and 3. We see good agreement with the theoret- 
ical predictions. In particular, data collapse is excellent 
(showing that the only memory left from the initial con- 




FIG. 3. tpi{t) (top) and Vtri{t) (bottom) versus t/r[m(0)]. 
To avoid overcrowding, only results for m(0) — —2/3, 0, and 
2/3 are shown. The dashed curves show the predictions 
l/2 + l/ln{t/T) and 0F/ln(t/r). 
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